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Abstract. I review the LHPC Collaboration's lattice QCD calculations of the generalized 
parton distributions of the nucleon and highlight those aspects of nucleon structure best 
illuminated by lattice QCD, the nucleon's spin decomposition and transverse quark structure. 



1. Introduction 

Generalized parton distributions provide the means to calculate aspects of nucleon structure not 
previously accessible to lattice QCD calculations. Primary among these are the nucleon spin 
decomposition 1 and the transverse quark structure The generalized parton distributions 
determine the fraction of nucleon spin carried by quark helicity, quark orbital motion, and 
gluons. Additionally, the generalized parton distributions provide a three dimensional picture of 
a fast moving nucleon by extending the ordinary parton distributions into the transverse plane. 

The Lattice Hadron Physics Collaboration is pursuing a program to calculate the nucleon's 
generalized parton distributions. Initial calculations have been completed for pion masses from 
750 MeV to 900 MeV [3 IH • Additionally current calculations are underway covering the region 
of pion masses from 350 MeV to 750 MeV In these proceedings I review the results from 
our earlier calculations highlighting the potential impact of lattice QCD calculations on our 
understanding of quark orbital motion and transverse quark structure within the nucleon. 

2. Generalized Form Factors 

Definitions of the generalized parton distributions can be found in 6 . However, there is 
an equivalent and more convenient language which is naturally suited to lattice calculations, 
generalized form factors. The generalized form factors are defined in terms of matrix elements 
of the twist two operators which can be calculated on the lattice. For example consider the 
tower of twist two operators O^ 1 " = qiD^ 1 ■ ■ ■ iD^™' 1 q. There is a corresponding tower 
of generalized form factors, A q ni , B^, and C^, defined by 

ri-l 

(P'\0^\P) = U(P')i E (Alm^ + B^K^+d^C^K^PiP) (1) 

i=0, even 

where the K are kinematic functions of P and P' and the spin labels have been suppressed. 



3. Lattice Calculation 

The details of our lattice calculation are given in [Jj and [2]. In brief, we calculate the matrix 
elements on the left hand side of equation ^ on the lattice. We then match these results to the 
MS scheme at fi = 2 GeV, and invert equation ^ to determine the generalized form factors. A 




key feature of our method is that we calculate an exhaustive set of matrix elements and solve 
the resulting overdetermined system of equations to extract the generalized form factors with 
as small an error as possible. In the following I use only results from a single lattice calculation 
with m n = 896(6) MeV. Further details can be found in the references above. Additionally, 
there are other similar calculations in the literature [SJ El 1101 111! I12| . 

4. Quark Orbital Motion 

The nucleon spin can be written as the sum of three separately gauge invariant observables, 

1/2 = 1/2 AY> u+d + L u+d + J 9 , (2) 

where the terms from left to right are the contributions from the quark helicity, quark orbital 
angular momentum, and total gluon angular momentum. The observables in equation |2] are 
each given by singlet matrix elements which pose a significant challenge to lattice calculations. 
Singlet matrix elements require the computation of disconnected quark diagrams which are 
computationally quite demanding. 

However non-singlet matrix elements receive no contributions from disconnected diagrams 
and are readily calculated on the lattice. In our calculation |13| . shown in figure Q we found 
that L u ~ d = —0.193(32) indicating that at least one of the quark flavors is undergoing significant 
orbital motion. This is quite interesting given that the connected contribution to L u+d is 0.002(3) 
which, ignoring disconnected diagrams, suggests a strong cancellation among the quark flavors. 
Due to the missing contributions to L u+d we are incapable of determining L u and L d separately, 
however, we can set a lower bound on the typical contribution from the quark orbital motion: 
the average magnitude of the quark orbital angular momentum, L avc = 1/2(|L"| + \L d \), is 
bounded from below, L ave > |L u_d |, indicating that on average the quarks are contributing, by 
way of their orbital motion, on the order of 20% of the nucleon's total spin. 

5. Transverse Quark Structure 

The transverse quark distribution, q(x,b±), gives the probability to find a quark q with a 
momentum fraction x at a distance 6j_ from the transverse center of the nucleon. The x moments 
of the transverse distribution are given by 

dx x"- 1 q(x, b ± ) = j !^e-^- A -Ay-Ai). (3) 



The Q 2 dependence of A^ is calculated directly on the lattice and, as explained shortly, for 
large n probes the transverse distribution of quarks with x near 1. Specifically we calculate the 
transverse rms radius and average x of the distribution / dx x n q(x, b±) allowing us to indirectly 
examine the x dependence of the transverse size of the nucleon. Additionally, upon entertaining 
a few assumptions, we can construct the b± dependence of q(x, b±) as well. The details of the 
following calculations are given in Jl] and [3]. 

5.1. Q 2 Dependence 

The Q 2 dependence of, and in particular the slope of, A q na determines the rms radius of the n th 
moment of q(x,b±), 

( b 2\ _ fd 2 b ± bif\dxx n - l q(x,b ± ) _ -4 dA^jO) 
\ ± '" Jd 2 b ± f\dxx^q(x,b ± ) A q n0 (0) dQ 2 

Notice that for large n the integrals in equation0]are dominated by the limit x — ► 1. Furthermore, 
as x — ► 1 a single quark carries all the longitudinal momentum of the nucleon. As discussed in [3], 
a quark in this limit is kinematically constrained to have b± — > 0. Additionally, explicit light 
cone wave functions show q(x,b±) behaves as S 2 (b±) for x — ► 1. Thus equation 0] demonstrates 
that the slope, and corresponding rms radius, of A^q vanishes asn-> oo. It is quite reassuring 
that we find, as shown in figure [2J that even for n = 1, 2, and 3 we see the slopes of A^q 
decreasing with increasing n indicating that the transverse distribution of quarks within the 
nucleon is indeed becoming more narrow as x approaches 1. 



5.2. x Dependence 

The rms radii {b 2 L ) n in equation 0] give the transverse size at a fixed moment n rather than at 
a fixed momentum fraction x, 

/ 2 v _ Jd 2 b ± 6j_g(x,6 ± ) 
\ fd 2 b ± q(x,b±) 

However for a fixed n there is a region of x which dominates Q?±) n , thus we can think of (&j_)„ 
as an average of over a region centered on the average x of the n th moment, 

_ Jd 2 b ± J\dx Ixlx^qjx,^) _ <x n > +2(-l) n Jd 2 b ± j^dxx n q(x) ^ < x n > 
n ~ fd 2 b ± J^dxx^qixA) ~ < x"- 1 > ~ <x n - 1 >' 

The anti-quark contribution in equation[5]is not accessible to current lattice calculations, however 
we expect it to be small and simple guesses based on the phenomenologically determined 
distributions indicated the following qualitative conclusion is not spoiled by ignoring it. In 
figure we plot y/(b 2 ± ) n versus (x) n demonstrating that the transverse size of the nucleon again 
appears to depend significantly on the momentum fraction at which it is probed. Note that 
(bj_) 2 i s n °t shown in the figure because n = 2 corresponds to the q + q distribution whereas 
n = 1, 3 correspond to the q — q distribution. 



5.3. b± Dependence 

Lattice calculations of parton distributions are restricted to a few low moments, however 
equation illustrates that the b± dependence of the lowest moments is quite illuminating. 
For n = 1 we have Jdx q(x,b±) which is the transverse probability distribution without regard 
to the momentum fraction. Furthermore the corresponding operator is conserved, hence q(b±) 




Figure 3. squares are (b\)^ d , triangles are 
{b\)^ d for n = 1,3 
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Figure 4. solid, dashed, dotted are q n (b±) 
for n = 1, 2, 3, q = u — d 



is independent of the renormalization scale. For n = 2 we have fdx x q(x,b±) which is the 
transverse momentum distribution for which the sum over all quarks and gluons is again 
conserved and hence the total transverse momentum distribution is also scale independent. 
Unfortunately the moments n > 3 no longer have simple physical interpretations, nonetheless 
they represent the transverse distribution of the corresponding moment. 

Assuming a dipole form for A^ , which is justified by figure El for Q 2 < 3.5GeV 2 , we can 
perform the integrals in equation |31 and calculate the b± dependence of each fixed moment for 
b± > (S.SGeV 2 )^ 1 . As figure illustrates we again observe that the lowest moment is the most 
broadly distributed in the transverse plane and higher moments are successively more narrow. 

6. Conclusion 

Our lattice calculations of the nucleon's generalized form factors have yielded insight into the 
quark structure of the nucleon. We have demonstrated that for a world with heavy pion masses 
the quark orbital motion constitutes roughly 20% of the nucleon's spin. Additionally we have 
shown, in several ways, that the transverse distribution of quarks within the nucleon depends 
significantly on the momentum fraction of the quarks. We are currently pursuing calculations 
with significantly lighter pion masses [U in an effort to examine the quark substructure of the 
physical nucleon. 
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